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Abstract. We prove a limit theorem for quantum stochastic differential equa- 
tions with unbounded coefficients which extends the Trotter-Kato theorem for 
contraction semigroups. From this theorem, general results on the conver- 
gence of approximations and singular perturbations are obtained. The results 
are illustrated in several examples of physical interest. 



1. Introduction 

It has been well established that the quantum stochastic differential equations of 
Hudson and Parthasarathy |18j provide an extremely rich source of realistic physi- 
cal models, particularly in the field of quantum optics 2 . The physical models in 
this framework are always Markovian, just like their classical counterparts, and it 
is the Markov property which makes these models particularly tractable. Several 
authors have developed mathematical methodology which can be used to obtain 
these Markovian models directly from fundamental models in quantum field theory 
by applying a suitable limiting procedure p] I14[ [6] , which places such models on 
a firm physical and mathematical foundation. On the other hand, the simplifica- 
tion provided by the Markov property is essential in many studies, particularly in 
investigations where nontrivial dynamics plays an important role. 

Despite the significant simplification provided by the Markov limit, many real- 
istic physical models remain rather complex, and often further simplifications are 
sought. In the physics literature it has been known for a long time that mod- 
els with multiple time scales can be significantly simplified, in parameter regimes 
where these time scales are well separated, by eliminating the fast variables from 
the description. Such singular perturbation methods are known as 'adiabatic elimi- 
nation' in the physics literature, and are extremely widespread on a heuristic level. 
Until recently, however, no rigoruous development was available in the quantum 
probability literature. In two recent papers [15, 3 J initial steps were made in this 
direction, but, as is discussed below, these results are not entirely satisfactory. 

The motivation for the current paper stems from the goal to develop a widely ap- 
plicable singular perturbation theory for quantum stochastic differential equations. 
As in previous results on the Markov limit, the simplified model is obtained by 
taking the limit of a sequence of quantum stochastic differential equations which 
depend on a parameter. Our main result, theorem [5] below, provides a general 
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method to prove the convergence of a sequence of quantum stochastic differential 
equations to a limiting equation. This technique hinges crucially on the Markov 
property of all the equations involved, and extends the Trotter-Kato theorem in the 
theory of one-parameter semigroups [5]. Using this technique, we develop a gen- 
eral singular perturbation method for quantum stochastic differential equations, 
theorem [3] below, which contains and extends all known results on this topic. 

Relation to previous work. To our knowledge, the first rigorous result on sin- 
gular perturbation for quantum stochastic differential equations was given in the 
paper [15] . This paper considers an atomic system inside an optical cavity, where 
the cavity is strongly coupled to an external field. It is shown that in the strong cou- 
pling limit, the cavity can be entirely eliminated from the model leaving an effective 
interaction between the atoms and the external field. The convergence to the limit- 
ing equation is given in the weak topology, and the proofs utilise Dyson expansion 
techniques similar to those used in the literature on Markov limits. However, the 
results in this paper are limited to the particular model under consideration, and it 
is unclear how to generalize these techniques to solve general singular perturbation 
problems. 

In [3j singular perturbations are studied in a general setting. The key observation 
in this paper is that, due to the Markovian nature of all the equations, the Trotter- 
Kato theorem can be used to prove strong convergence. This is done by defining 
certain skew semigroups on the algebra of bounded operators on the initial space. 
Though results of a general nature are obtained in this fashion, the technique is 
essentially restricted to the case that all quantum stochastic differential equations 
of interest have bounded coefficients. In the unbounded case the semigroups on 
the algebra of bounded operators are typically only weak* continuous and not 
strongly continuous, which precludes the application of the Trotter-Kato theorem. 
Unfortunately, the boundedness assumption excludes a large number of singular 
perturbation problems of practical interest, including the problem studied in [T5] . 

The singular perturbation results in this paper extend and bridge the gap be- 
tween the results of [H2 [3] . We consider a family of quantum stochastic differential 
equations whose coefficients are not necessarily bounded but possess a common 
invariant domain. Rather than considering semigroups on the algebra of bounded 
operators, we work directly with semigroups on the initial Hilbert space. This 
approach has many far reaching consequences. Unbounded coefficients are easily 
treated in this setting, and the type of convergence which we obtain (strong conver- 
gence uniformly on compact time intervals) is much stronger than in [151 13] ■ More- 
over, the proofs in this setting are significantly simplified as compared to |15[ 15]. 
which highlights that the current approach is very natural. In order to prove our 
singular perturbation results, we first prove a general theorem on the convergence 
of quantum stochastic differential equations with unbounded coefficients, which ex- 
tends the Trotter-Kato theorem to our setting. With this theorem in hand, the 
remaining work on singular perturbations is chiefly algebraic (rather than analytic) 
in nature. 

Though motivated by a rather different set of problems, our main convergence 
result, theorem [2j is closely related to results obtained for the purpose of proving 
existence and uniqueness of solutions of quantum stochastic differential equations 
with unbounded coefficients [HJ proposition 3.4] and [211 theorem 1.7], [20]. Instead 
of showing the existence of a contractive cocycle satisfying the limit equation, we 
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assume the existence of a unitary cocycle satisfying the limit equation, which allows 
us to obtain a much stronger form of convergence than is obtained in [ITJ EH EQ] . 



Organization of the paper. The remainder of this paper is organized as follows. 
In section[5]we introduce the class of models that we will consider and the technical 
conditions which are assumed to be in force throughout the paper. In subsection l2.ll 
we describe our main convergence result, and subsection 12.21 describes our general 
result on singular perturbations. Section [3] is devoted to the discussion of several 
examples of physical interest which demonstrate the flexibility of our results. Fi- 
nally, section 0] presents the proofs of the results described in subsection 12. 11 while 
section [5] presents the proofs of the results described in section 12.21 

2. Main results 

Throughout this paper we let 7i, the initial space, be a separable (complex) 
Hilbert space. We denote by T — r s (L 2 (R + ; C™)) the symmetric Fock space with 
multiplicity nfN (i.e., the one-particle space is C™ <g> L 2 (R + ) = L' 2 (R + ; C™)), and 
by e(/), / € L 2 (R + ;C n ) the exponential vectors in T. The annihilation, creation 
and gauge processes on T, as well as their ampliations to Ti <g> J 7 , will be denoted 
as A\, A 1 ^ and K\\ respectively (the channel indices are relative to the canonical 
basis of C"). Moreover, we will fix once and for all a dense domain D C 7i and 
a dense domain of exponential vectors £ = span{e(/) : / S &} C where 
6 C L 2 (K + ; C n ) R L^ C (M.+ ; C") is an admissible subspace in the sense of Hudson- 
Parthasarathy [18 which is presumed to contain at least all simple functions. For a 
detailed description of these definitions and of the Hudson-Parthasarathy stochastic 
calculus which we will use throughout this paper, we refer to [TB I 122" ! 12], 

For every k € N we consider a quantum stochastic differential equation of the 
form 
(1) 

!n n n 

£ - dA ? + £ M i k)dA ? + £ 4 k)dA i + K(k)dt 
i,j=l i=l z=l 

(k) 

where U Q = I and the quantum stochastic integrals are defined relative to the 
domain T>®£. The purpose of this paper is to prove that, when the dependence of 
the coefficients on k is chosen appropriately, the solutions converge as k — * oo 
in a suitable sense to the solution of a limit quantum stochastic differential equation 

{n n n 

£ (Ny - ) dAl 3 + Mi dA? + ^ L, dA\ + K dt 
— l i—1 i—1 

Of course, the form of K, etc. will depend on our choice for K^ k \ etc., and a part 
of our task will be to identify the appropriate limit coefficients. We will both prove 
a general convergence result — a version of the Trotter-Kato theorem for quantum 
stochastic differential equations — and investigate in detail a large class of singular 
perturbation problems (known as adiabatic elimination in the physics literature) 
which are of significant practical importance in obtaining tractable simplifications 
of complex physical models. 

(k) 

As it turns out, it is not always the case that J7 t converges on the entire Hilbert 
space Hd^F; in singular perturbation problems the limit will only exist on a closed 
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subspace Ho <8> T, while for vectors in Hq ® T the limit becomes increasingly 
singular as k — > oo. We will return to this phenomenon in detail later on; for the 
time being, suffice it to say that we must take into account the possibility that the 
limiting equation ([2]) lives on a smaller space Ho ® -F than the prelimit equations 
(JTJ. For the time being, we will fix the closed subspace Ho C H and the dense 
domain T>q C T> fl 7io in Wo- The quantum stochastic integrals in ([2|) will then be 
defined relative to the domain T>q®£. 

We begin by imposing conditions on these equations that are assumed to be in 
force throughout this paper. For notational simplicity, we use the same notation 
for operators on H or on J- and for their ampliations to H ® T (and similarly for 
Ho and Ho ® J 7 , etc.) Throughout, we denote by t an adjoint relationship on the 
domain of definition, i.e., two operators X, X' on the domain D (or 2? ) are always 
presumed to satisfy the adjoint pair property (u,Xv) = {X'u,v) for all u,v £ T> 
(resp. re- 
condition 1 (Hudson-Parthasarathy conditions). For all k £ N and 1 < i,j < n, 
the operators K ik \ if( fc )t, hf\ L (k) \ M^ ] , , , have the com- 

mon invariant domain D, and the operators K, K\ Li, L\, Mi, M\, N^, nJj 
have the common invariant domain T>o- In addition, we assume that the Hudson- 
Parthasarathy conditions 

n n 



Y,N mj Nt j =Y J N} m N j t = 5 ml 

3=1 3 = 1 

are satisfied on T>o, and that 

n n 



3=1 



3 m 

3=1 3=1 

are satisfied on T> for all k £N. 

Denote by 6 t : L 2 ([t, oo[; C") -> i 2 (R+; C") the canonical shift 9 t f(s) = f(t+s), 
and by 0t : J-[ t — > .F its second quantization (here T = JF t j ® jFj t denotes the 
usual continuous tensor product decomposition). Recall that an adapted process 
{Ut ■ t > 0} on H ® J 7 is called a contraction cocycle if C/j is a contraction for all 
t > 0, < JJ t is strongly continuous and /7 s+t = C/ s (/®9*C/ t 6 s ) with 7(g)9*C/ t 0s € 
-F s ] ®H®T[ S = H®T. If C/ t , i > are even unitary, the process is called a unitary 
cocycle. 

Condition 2 (Cocycle solutions). For a/Z k £ N, equation {!]) possesses a unique 
solution {f/ t : i > 0} which extends to a contraction cocycle on H ® T ' , while 
equation (0) possesses a unique solution {Ut '■ t > 0} which extends to a unitary 
cocycle on Ho ® -F. 

The following elementary result is proved in section 2] 
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Lemma 1. For a, f3 € C n , define T t : Hq — > 7io such that 

(u,T t {al3) v) =e- (|a|2+l/3|2)t/2 ( M ®e(«/ [0 , t] ),C/ tW ®e(/3/ [0 , t] )) V«,i;e7io, t > 0. 



TTien T t is a strongly continuous contraction semigroup on Ho, and the generator 



The same result holds for j 1 ^' : H — > H and j£( k > a P) f defined by replacing Ut by 
and making the obvious modifications. In particular, Dom(«Sf z) T>. 

We now impose one further condition: we assume that §5§ completely determines 

_£>(a/3)_ 

Condition 3 (Core). For any a, ft G C™, the domain T>o is a core 

(k) 

Remark 1. Note that it is not necessary for our purposes to assume that t/ t 
is unitary or that I? is a core for JfC^/ 3 ), fceN. Typically this will nonetheless 
be the case in physical applications; in particular, recall that when the coefficients 
are bounded, condition Q] guarantees that the solution of the associated quantum 
stochastic differential equation is unitary and the core condition is trivially satisfied. 
In the unbounded case, however, the solution may only be a contraction due to the 
possibility of explosion. 

Remark 2. We have chosen the right Hudson-Parthasarathy equations |T]) and 
(J2J), rather than the more familiar left equations where the coefficients are placed 
to the left of the solution. This means that the Schrodinger evolution of a state 
vector tp 6 Ho (8 T is given by U^if, etc. The main reason for this choice is 
that for quantum stochastic differential equations with unbounded coefficients, it 
is generally much easier to prove the existence of a unique cocycle solution for the 
right equation than for the left equation (see [13110]); this is chiefly because it is not 
clear that the solution should leave T>o®£ invariant, so that the left equation may 
not be well defined (this is not a problem for the right equation, as the solution only 
appears to the left of all unbounded coefficients and is presumed to be bounded). 
This appears to be an artefact of the fact that quantum stochastic integrals are 
defined on a fixed domain, and is not necessarily a physical problem. 

We will ultimately prove that C/ t converges to C/ t * strongly on Ho® J 7 uniformly 
on compact time intervals, i.e., we will show that 



Therefore, there is no loss of generality in working with the more tractable right 
equations. If we wish to begin with a well defined left equation, our results can be 
immediately applied to the Hudson-Parthasarathy equation for its adjoint. 

Remark 3. In practice, the result of Fagnola [10l [12] provides a convenient suf- 
ficient condition which can be used to verify conditions [2] and [3] Let us recall a 
slightly stronger version of this result (for simplicity phrased in terms of Ut, the 
analogous result holds for C/ t ). Suppose that for all u € T>o and I E N, there 
exists a constant c(£, u) such that: 




k^oo o<t<T 



lim sup \\Uf'*i/) - U?ip\\ = 



Vi/iGWo®^ T < oo. 
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(a) For all u £ T>q and for some e > independent of u 

oo 

^^c{l,u)e P ' < oo; 
1=1 

(b) For all I £ N and all choices of X(l), . . . , X(£), where each X(-) is one of 
K, K\ L u L\, M u Mj, Nij, N^, we have for all u £ V 

\\X(1) ■ ■ ■ X(£)u\\ < c(£, u) y/(e + m)\, 

where m is the number of occurences of K or K' in the sequence X (1) , . . . , X(£) . 

Then equation J2J) possesses a unique solution {Ut ■ t > 0} that extends to a unitary 
cocycle on Ho0!F, which verifies condition [5] To show that T>q is a core 
it suffices to note that the Fagnola conditions imply that every u £ Vq is analytic [H 
definition 3.1.17], so that the result follows from [4j corollary 3.1.20]. This verifies 
condition [3] 

If we are only interested in the existence and uniqueness of a contraction cocycle 
Ut, simple sufficient conditions exist in a much more general setting than the one 
considered in lO j [T2 | . In particular, from [201 theorem 1.1], one can see that it 
suffices to prove that the closure of Jz? on Vq, as defined in equation ((3|), is the 
generator of a strongly continuous contraction semigroup on TLq for every a, £ C™. 

2.1. A TVotter-Kato theorem for quantum stochastic differential equa- 
tions. In order to prove that converges to Ut as k — > oo, we would like to 
use an argument which requires only to verify conditions on the coefficients of the 
associated equations. After all, the dependence of the coefficients on k is known 
in advance, while the solutions of the equations could potentially depend on k in 
a complicated manner. In order to set the stage for our main convergence result, 
let us shift our attention for the moment to the semigroups T^ k ' a ^ and T^ a ^ . 
The convergence of the former to the latter can be verified using a standard tech- 
nique: the Trotter-Kato theorem. For future reference, we state this result here in 
a particularly convenient form [SI theorem 3.17, p. 80]. 

Theorem 1 (Trotter-Kato). Let Ti be a Hilbert space and let H.q C TL be a closed 

(k) 

subspace. For each k £ N, let T t be a strongly continuous contraction semigroup 
on Ti. with generator JfW . Moreover, let T t be a strongly continuous contraction 
semigroup on H.q with generator _Sf . Let T>q be a core for Jzf. The following condi- 
tions are eguivalent: 

(a) For all if} £T> there exist ip^ £ Dom(Jf(' c )) such that 

(b) For all T < oo and tp £ Ho 

lim sup \\T t (k) ijj - T t ip\\ = 0. 

k^oo o<t<T 

It should be noted, in particular, that the condition (a) depends directly on the 
generators and J>? , whose form (in terms of k) is typically known explicitly. 

Condition (b), however, entails the convergence of the semigroups (in a particularly 
strong form). 
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The central result of this paper is an extension of the Trotter-Kato theorem to 
the quantum stochastic differential equations (QSDEs) defined in equations ^ and 
@. The proof of the following theorem can be found in section 0J 

Theorem 2 (QSDE Trotter-Kato). The following conditions are equivalent. 

(a) For every a, (3 € C™ and u 6 T>q, there exist vS k > 6 Dom(^f W'W) such that 

(b) For every tp € Hq <8> J 7 fl*"^ T < oo 

lim sup ||[/ t (fe) > - = 0. 

Moreover, ifD is a core for all^ af3 \ k € N i/iera we caw always choose {u' fe '} C 
P. 

The convergence of T^ k ' a ^ to T^ a ^ for all a, /? £ C n thus turns out to be equiv- 
alent to the convergence, in a very strong sense, of to C/ t . This is a powerful 
technique precisely because of the fact that convergence can be determined directly 
from the generators ^f( k:a ^ , jgf ( Q ^) which depend only on the coefficients of equa- 
tions (Q]) and ||2J); the technique does not require us to understand precisely the 

(k) 

way in which U{ depends on k. The benefit of this approach is clearly demon- 
strated by the proof of the singular perturbation results below, which follows as an 
application of theorem [2j 

Remark 4. Note that for the implication (a)=>(b) it always suffices to choose 
{u^} C 2?, as X> C Dom(:Sf (fc;a/J )). It is only when we wish to ensure the existence 
of {u^} C V in the reverse implication that we must assume T> to be a core for 
all &( k > a P) , k £ N. In practice we are almost always interested in (a)=>(b), where 
this is irrelevant. 

Remark 5. The proof of theorem [2] is easily modified to show that \\uj: i[> — 
UtipW k ^°° ) for every ip £ TIq ® T and t < oo. However, uniform convergence on 
compact time intervals appears to be easier to obtain for convergence to the adjoint 
U£. As it is ?7 t * which defines the physical time evolution, we do not attempt to 
obtain uniform convergence to Ut- 

2.2. Singular perturbations. Let us start by introducing the basic assumptions. 
We consider Ti and D as fixed from the outset; as we will shortly see, we must 
choose 7io CWina particular way. 

Assumption 1 (Singular scaling). There exist operators Y, Y>, A, A\ B, B* , Fi, 

F},Gi, Ot, Wij, W} 3 with the common invariant domain T> such that 

K (k) = k 2y + kA + Bi L (k) = m + N (k) = w .. 

for all k € N and 1 < i, j < n. 

Remark 6. Note that condition [T] imposes additional assumptions which are pre- 
sumed to be in force throughout. In particular, the Hudson-Parthasarathy relations 
in condition [I] determine completely the form of the fourth coefficient . 
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The interpretation of this scaling is that the dynamics of the physical system 
contains components that evolve on a fast time scale and on a slow time scale. The 
dynamics on the fast time scale is generated by quadratic (oc k 2 ) term, while the 
coupling between the fast and slow time scales is determined by the linear (oc k) 
term. As k — > oo, the dynamical evolution of the fast components of the system 
will thus become increasingly singular. However, under appropriate conditions, the 
dynamical evolution of the slow components of the system will converge as k — > oo, 
and the time evolution of these components in the limit will be decoupled from 
the fast components. Thus, in the limit of infinite separation of time scales, we 
have 'adiabatically eliminated' the fast components of the system leaving only an 
effective time evolution of the slow components. Typically the slow components of 
the system are the physical quantities of interest, and the elimination of the fast 
components from the model constitutes a significant simplification in many complex 
physical models. 

As only slow dynamics should remain in the limit, we aim to choose Ho such that 
it contains only the slow degrees of freedom of the system. In order to ensure that 
the slow dynamics do indeed have a well defined limit, we must impose suitable 
structural assumptions on the coefficients; in particular, we must enforce that the 
terms which are linear and quadratic in k do not directly generate dynamics on 
Ho — in the latter case, the limit would surely be singular on Ho- We presently 
collect all the necessary structural assumptions; their significance can be clearly 
seen in the proof of our main result. 

Assumption 2 (Structural requirements). There is a closed subspace Ho C H 
such that 

(a) P V C V; 

(b) YP = on V; 

(c) There exist Y , Y' with the common invariant domain T>, so that YY = 
YY = P l ; 

(d) F]Pq = onV for all l<j< n; 

(e) PqAPo = on V. 

Here Pq denotes the orthogonal projection onto Ho and P\ the orthogonal projection 
onto H.Q . We choose the dense domain T>o = PqD in Ho- 

It remains to introduce the limit coefficients. The following expressions emerge 
naturally in the proof of our main result through the application of theorem [2j 

Assumption 3 (Limit coefficients) . Define the operators on Ho 

K = P (B - AYA)P Q , U = Po(G t - AYFi)P , 

n n 

M i = ~Y J PoWij(G] - f]yA)P , N tj = PoWuiFjYFj + %)P . 
j=i t=i 

Then K , , Li, h\, Mi, Mj , Nij, Nj- have common invariant domain T>o- To en- 
sure that these coefficients satisfy the Hudson-Parthasarathy relations of condition 
[H we require 

n n 

P (G t - AYFi)Pi =J2 p oWu(FjYF j + %)P X = £ P 1 W u (F^YF j + %)P = 0. 
e=i i=\ 
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Our first order of business is to verify that the coefficients K, Li, Mi, Nij do 
indeed satisfy the Hudson-Parthasarathy relations in condition [TJ This is proved in 
section [5J 

Lemma 2. Under assumptions QH3 the coefficients K, Li, Mi, Nij satisfy the 
Hudson-Parthasarathy relations in condition^ 

Finally, let us state our main result, whose proof will also be given in section [5] 
We remind the reader that beside the above assumptions, conditions [T]-[3] are still 
presumed to be in force; in particular, one must verify separately that ([T]) and @ 
uniquely define contraction and unitary cocycles and that Dq is a core for J2f( a ^), 
a, [3 e C n . 

Theorem 3 (Singular perturbation). Under assumptions QH3 the singularly per- 
turbed equations (QP converge to the limiting equation JJ|) on H.q: 

lim sup \\U$ k) *ip-UZil>\\=Q for all tji € H ® J 7 . 

fc^oo o<t<T 

3. Examples 

3.1. Atom-cavity models with bounded coupling operators. In this sub- 
section, we will consider a class of physical models which consist of a harmonic 
oscillator coupled to auxiliary degrees of freedom. Such models cover various ap- 
plications in quantum optics, such as a single mode optical cavity coupled to the 
internal degrees of freedom of a collection of atoms. Set Ti = Ti' <g> £ 2 (Z + ). On 
the canonical basis {ipi : i € Z + } of ^ 2 (Z + ), the creation, annihilation and number 
operators are defined as 

ipi. = \Ji + 1 tpi+i, bipi = y/iipi-i, Wbipi=iipi. 

These definitions extend directly to the domain V = span{</?i : i 6 Z + } C £ 2 (Z + ), 
and we will choose the dense domain T> = TC ® T>' 'mTt. We now consider prelimit 
equations |lj whose coefficients have the following form on D: 

N f = , L<*> = F^ fot +G f ) ; K W = E W b i b+E [V h \ +E W b + E W . 

Here Epq , F^ k \ G\ k \ are bounded operators on Ti' , and condition Q] is pre- 

sumed to be in force (in particular, M\ are completely determined by these defi- 
nitions) . 

We begin by showing that such equations possess unique solutions which extend 
to unitary cocycles on Ti ® T, and even that I? is a core for jjfC ;"/ 3 ) , a,/3 € C. 
We will establish these facts by verifying the conditions of Fagnola, see remark [3] 
Indeed, this is only a minor modification of the computations performed in [101 112j , 
and we refer to those papers for a more detailed account of the necessary steps. 

Lemma 3. The conditions of Fagnola for the existence of unique unitary cocycle 
solutions (remark^) are satisfied for the class of models considered in this subsec- 
tion. 

Proof. Note that all coefficients are linear combinations of operators of the form 
X¥b, Xb\ Xb, X, where X are bounded operators on Ti' whose norms are 
bounded by a fixed constant ||X|| < K < oo (for simplicity, choose K > 1). 
Arguing as in [101 02] , we find 

\\X(l)---X(e)iP®<p p \\ <pV(p + * + m)!, 
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where tp <E H', X(l), . . . X(£) is an arbitrary selection of coefficients or their ad- 
joints, to is the number of occurences of K ( fe ) or ifWt, and p = AK 2 . Using the 
elementary estimate (£ + m)\ < 2 e+m £\m\, we can thus choose (see remark [3] for 
notations) 

c(4V®¥> P ) = p e 2 e ^/pJ¥, 
which does indeed satisfy the necessary analyticity requirement for sufficiently small 
e. But any vector u S T> is the linear combination of a finite number of vectors of 
the form ip ® ip p , so that the claim follows easily. □ 

With the existence and uniqueness taken care of, it becomes straightforward to 
apply our main result on singular perturbations, theorem [3l We will demonstrate 
this procedure in two physical applications: the adiabatic elimination of a single 
mode cavity, which extends the results of Gough and Van Handel [15], and the 
simultaneous elimination of a single mode cavity and the excited level of an atom 
as considered in Duan and Kimble [8j. 

Example 1 (Elimination of the cavity). We consider the case where the coupling 
of the oscillator with the external field becomes increasingly strong; in the limit, we 
expect to eliminate the oscillator entirely from the model as it will be forced into 
its ground state. In quantum optics, this corresponds to the 'adiabatic elimination' 
of an optical cavity in the strong damping limit. To make this precise, we set 

ij ~ i ~ *> i ~~ *' 

E^l = k 2 Eu, — kE w , E^ = kEoi, E$ = E o 

for 1 < i,j i < n, where where Sij, Fi, Gi, E pq are bounded operators on 7i which 
are chosen such that the Hudson-Parthasarathy relation of condition [1] are satisfied 
for all k E Af. We obtain the following result, which does not presume conditions 
[U and El 

Proposition 1. Suppose that En has a bounded inverse. SetTio = T>o = H'^dpo, 
K = E ao — E i(En) 1 E W , Li — d — Eqi(Eh) 1 F i , 

n 

N ij = Y / S i e{FnE n )- 1 F j + 6z j } 

on T>q , and define Mi through the Hudson-Parthasarathy condition [3 Then 
lim sup \\U[ k) *i) - U*i>\\ = for all i> £ H O T. 

fc^oo o< t <T 

Proof. Note that condition [T] is satisfied by assumption for the prelimit equations 
(IlJ, while condition[2]is satisfied for the prelimit equations by the Fagnola conditions 
verified above. Note that the singular scaling assumption [1] is also satisfied by 
construction. Next, we turn to the structural requirements of assumption O It is 
immediate that Vq = PqD c £>, while Y = En b^b evidently vanishes on T>q. We 
define Y on T> by setting Y tp® (po — and Y ip® p p = p^ 1 (En)~ 1 ip^ ip p for p > 1; 
then indeed YY = YY = Pi on V. It is immediate that Fjb vanishes on Vq and that 
A = E\o w + Eq\ b satisfies PqAPq = 0. Hence assumption [2] is satisfied. Next we 
note that, taking into account the identity X'bYb^X ip (E> <fio = X' (En)' 1 X ip ® (po 
for any bounded operators X, X' on H' , the coefficients defined in the proposition 
coincide precisely with the coefficients in assumption 02 and it is easily verified that 
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the remaining conditions of assumption [3] are also satisfied. In particular, the limit 
coefficients satisfy the requirements of condition [TJ But as the limit coefficients are 
bounded, conditions [5] and [3] are automatically satisfied. We have thus verified all 
the requirements of theorem [3l □ 

Example 2 (Duan-Kimble). We consider a three level atom coupled to a cavity 
which is itself coupled to a single output field, i.e., we set Ti' = C 3 and n = 1. Let 
us denote by (|e), |+), |— )) the canonical basis in C 3 . In this basis we define 

/0 I 0\ /0 1 N 

cr^ +) = , 0-^=000 
\0 0/ \0 y 

Define moreover = a^* and P± = The quantum stochastic differ- 

ential equation for a lambda system with one leg resonantly coupled to the cavity, 
under the rotating wave approximation and in the rotating frame, has the following 
coefficients: 

= -Ef> = k 2 ga { +\ E$ = k(a^a* - a^a), 

where 7 £ M. + and a S C. This model coincides with that of Duan and Kimble 
[5] , except that we have allowed for additional driving on the uncoupled leg of the 
atom with amplitude a. Let us now define operators Y, A, B, F, G, W as 

Y = — — fetft + g (o-Wftt _ .Wft) j A = a^a* - a^a, 

B = 0, F 1 = ^b\ G 1= 0, W n =I. 

These definitions are easily verified to satisfy condition [T] and assumption [TJ 

We proceed to verify assumption [2j Set Ho = V = span{|+) <g) ip , \— ) ® ip }- 
Then it is easily verified that V = P V C V and that yP = F^Pq = P AP = 
on T>. It remains to define Y. To this end, define the following subspaces of Ti: 

Hj = span{|+) <gs (fj, \-)<8><fj, \e)®ipj-i}, j € N. 

Note that Pi I? = and that the subspaces Hj, jeN are all invariant under 

the action of Y where, with respect to the basis (|+) <£> (fj, |— ) ® <^j, |e) ® <Pj-i), 

We may then construct Y = ©Jl 1 (y|w j ) _1 on P{D (and yX>o = 0), which reads 




' 7(3-1) 



<?v7 A 



^=4 o^o, rf ,=i4^i + , 2j 



TJ 
2 

in the same basis as above. Thus assumption [2] has been verified. We are now in 
a position to introduce the limit coefficients. By assumption [3J we must evidently 
define on Tio 

K = -^P-, L 1 = -^i->4+\ N U =I-2P_. 

V 9 
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It is again easily verified by explicit computation that the remaining conditions of 
assumption [3] are satisfied. Moreover, all the coefficients that we have introduced 
satisfy the Hudson-Parthasarathy relations of condition [TJ condition [5] is satisfied 
for the prelimit equations by the Fagnola conditions, and conditions [T] and [3] are 
automatically satisfied as the limit coefficients are bounded. Thus all the conditions 
of theorem [3] have been verified. 

3.2. Coupled oscillators. In the examples in the previous subsection, the appli- 
cation of theorem [3] was significantly simplified by two convenient properties: the 
Fagnola conditions could be verified for that class of models, and the limit equa- 
tions always had bounded coefficients. In the present section, we will develop an 
example which enjoys neither of these properties. 

We consider a single mode cavity in which one of the mirrors is allowed to 
oscillate along the cavity axis (see, e.g., [7]). The initial system will therefore have 
two degrees of freedom: the kinematic observables of the oscillating mirror and 
the usual observables associated with the cavity mode. This is implemented by 
choosing H = L 2 (M) <g) £ 2 (Z + ) . For the cavity mode Hilbert space £ 2 (Z+ ) , we define 
the domain V of finite particle vectors and the creation and annihilation operators 
b\ b as in the previous section. In the mirror Hilbert space L 2 (R) we will choose 
as dense domain the Schwartz space [13l ch. 8-9] 

S = {/ e C°°(R) : WD^flU < oo Va,/3 G Z+}, 

where we have defined the operators 

Df(x) = $P-, Xf(x)=xf(x) V/eC°°(R). 
ax 

Note that these operators leave S invariant, as do all operators of the following 
form: 

(g(X)f)(x) = g(x)f{x) VfeS,ge C°°(R) s.t. Halloo < oo for all (3 e Z+. 

In addition, let us recall the following well known facts [23l appendix V.3]. If we 
define the operators B, B^ with invariant domain S and the vectors $j £ S as 

b = st = = n-^ur^ye-^ 2 , j e z + , 

then S' = span{$j : j e Z + } C S is dense in i 2 (M), B^BQj = B$j = 

\fj i, and B'&j = y/j + 1 Therefore S' corresponds to the domain of fi- 

nite particle vectors in ^ 2 (Z + ) in the usual isomorphism between L 2 (M) and £ 2 (Z + ). 
In the following we will set T> = S ® V, but we will occasionally exploit S' whenever 
this is convenient. 

We are now in a position to introduce the prelimit equations |T]) . We set n = 1, 
i.e., there is only one external field, and choose the following prelimit coefficients: 

N$=I, L[ k) =k^b\ K (k) =k 2 jtf(B + £t)_2 tfb + i,VLB ] B, 

where 7, Q > are fixed parameters. The physical interpretation of this model 
is that the optical frequency of the cavity mode is determined by the length of the 
cavity, and hence by the displacement (B + B') of the mirror; this accounts for the 
(B + B*)tfb term. At the same time, we presume that the mirror is oscillating in a 
quadratic potential, which accounts for the B^ B term, while the remaining terms 
couple the cavity mode to the external field. In the strong damping limit k — > 00 
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we expect as before that the cavity is eliminated, so that we obtain an effective 
interaction between the mirror and the external field. 

Lemma 4. The prelimit equations (QP each possess a unique solution which extends 
to a contraction cocycle onTi® T ' . 

Proof. First, note that condition [T] is satisfied for these equations on the domain T>. 
It suffices to verify that the closure of Jz? ^'• a P s > on T> is the generator of a strongly 
continuous contraction semigroup for every k £ N and a, (3 £ C; see remark [3] and 
[20] theorem 1.1]. Fix k,a,f3, and denote by Jz? the operator ££( k ;°-P) with the 
domain D. By an easy corollary of the Lumer-Phillips theorem [20J lemma 2.1], it 
suffices to prove that D is a core for j£f*. But this follows immediately from the 
following lemma, and the proof is complete. □ 

Lemma 5. Let T be an operator on the domain T> of the following form: 

T = ci + c 2 6 f 6 + c 3 B^B + c 4 b + c 5 6 f + c 6 {B + S t )6 t &, 
where Ci g € C. Then T> is a core for T* . 

Proof. As (u,Tv) = (T^u, v) for u,v £ V, clearly V C Dom(T*) where T*\ v = Tt. 
Our goal is to show that the closure clT> of T' is in fact T*, i.e., we must show 
that Dom(T*) C Dom(clT t ). We begin by noting that 



OO 

Dom(clT t ) = \ ^ dj,e $j O w : ^ |<2j,.e| 2 < oo and ^ \Tj,e(d)\ 2 < oo 
j,i=o j,e=o j,e=o 



where we have written 

Tj, t {d) = (4 + c* 2 e + 1,/,., + ctVidj,i-i 



+ 4Vi + 1 d Jlt+ i + cg(v7+ldj+M + V7 d,-M)£ 

Indeed, this space is evidently the completion of S' ® 2?' with respect to the Sobolev 
norm ||i>||Tt = IMI + \\T *v\\, and V is already included in this space (this follows 
from the representation [23, theorem V.13]). It thus remains to show that Dom(T*) 
is included in this space as well. To this end, let us define for all N £ N 



^0 ^0 IIE^=0^W^®Wll 

Thus || tt at || = 1 and um £ T) for all JV e N. Now recall that by the Riesz rep- 
resentation theorem v £ Dom(T*) iff \{Tu,v)\ < C\\u\\ for all u £ V. But if 
v <£ Dom(clTt), then 



\{Tu N ,V)\ = \(Tu N ,P N+1 v)\ = \{u N ,T r P N+1 v}\ 



N 
3,1=0 



TV— >oo 
► OO 



where Pn is the orthogonal projection onto the subspace span{<I>j <g) (pg : j,l < N}. 
Thus v Dom(clTt) implies v Dom(T*), so that Dom(T*) C Dom(clT+) as 
desired. □ 

Let us now turn to the assumptions HHHl Clearly our coefficients are of the form 
required by assumption!]] Let us set Ho = L 2 (R) <g> Ct^o and T> = P V = S®Ctpo. 
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Then all the requirements of assumption [2] are clearly satisfied except that we must 
verify that 

Y = id(B + B^) - ~ tfb 
possesses an inverse Y on P{D. We may simply set, however, 

1 r + 7i 

Y ip <g> (po = 0, Yip®tp j =j 1 \id{B + B^) - ij)®<pj, j>l. 

As B + B^ = XV2 and the function g : K -> C, = (ifixVZ - 7/2)" 1 is smooth 
and bounded together with all its derivatives, the bounded operator Y leaves the 
Schwartz space invariant, and clearly YY = YY = Pi on T>. Therefore, all the 
conditions of assumption [2] arc satisfied. It should be noted that it would not have 
been sufficient to choose the smaller finite excitation domain S' in L 2 (R), despite 
the fact that it is left invariant by all prelimit coefficients, as the inverse Y docs 
not leave S' invariant. 

Having verified assumptions Q] and [21 the form of the limit equation @ is deter- 
mined by assumption [3J In particular, we must choose the limit coefficients 

M? OB + £t) + 7/2 + 

N " = mB + Bi)-]/2 > Li = ' K=inBiB 

on T>o, and one can verify by straightforward manipulation that the remaining 
conditions in assumption [3] are satisfied. However, we have not yet established that 
the limit equation possesses a unique solution which extends to a unitary cocycle, 
and whether the core condition [3] is satisfied. This is our next order of business. 

Lemma 6. The limit equation (0) possesses a unique solution which extends to a 
contraction cocycle on Tl ® T . Moreover, the core condition^ is satisfied. 

Proof. First, note that condition [1] is satisfied for the limit equation on the domain 
V . It suffices to verify that the closure of Jzf on T> is the generator of a strongly 
continuous contraction semigroup for every a, f3 € C: existence and uniqueness 
follows from remark [3] and [5D1 theorem 1.1], while the fact that 2?o is a core for the 
corresponding generators is then trivially satisfied. Note that we have 

^ = (x* j3 + +7/2 _ H!M + m B^B. 

H i${B + St) - 7 /2 2 

But as B^B is essentially self-adjoint on S, the closure of zO B^B on S is the gener- 
ator of a strongly continuous unitary group (and hence of a contraction semigroup) , 
while the first two terms in the expression for ( a ^) constitute a bounded and dis- 
sipative perturbation. The result therefore follows from a standard perturbation 
argument [5l theorem 3.7]. □ 

Lemma 7. The solution of equation fl|] extends to a unitary cocycle. 

Proof. As we have shown that there is a unique solution, it suffices to show that 
there exists a unitary process Ut that satisfies equation @. To this end, consider 
the equations 

S t = I+ f SsiilB^Bds, R t =I+ f R. (S.NuSt - I) dA„. 
Jo Jo 

Then St is the strongly continuous unitary group generated by the closure of iflB^B, 
while R t satisfies a Hudson-Parthasarathy equation with bounded, though time 
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dependent, coefficients. The latter possesses a unitary solution as can be verified, 
e.g., using the results in [17]. Moreover, using the quantum stochastic calculus, it 
is immediately verified that Ut — RtSt satisfies equation and RtSt is clearly 
unitary. □ 

We have finally verified all the conditions HHS] and assumptions [THIS We can thus 
invoke theorem G2 and we find that indeed 

lim sup \\U [ t k) *i> - U*Tp\\ = for all tp £ H <8> T 

fe^oo o<t<T 

with the prelimit and limit coefficients as define above. 

3.3. Finite dimensional approximations. Suppose that we are given a quan- 
tum stochastic differential equation of the form @ with dim7io — 00 ■ Though such 
an equation may be a realistic physical model, it can not be simulated directly on 
a computer. To perform numerical computations (typically one would simulate a 
derivative of this equation, such as a master equation or a quantum filtering equa- 
tion), we must first approximate this infinite dimensional equation by one which is 
finite dimensional. A very common way of doing this is to fix an orthonormal basis 
{^i}e>o C 2?0i so that we can approximate the coefficients in the equation for Ut by 
their truncations with respect to the first k basis elements. We will show that the 
solutions of the truncated equations do in fact converge to Ut as k — > oo. Though 
the result is of some interest in itself, it also serves as an exceedingly simple demon- 
stration of the Trotter-Kato theorem [2] which is not of the singular perturbation 
type (theorem [3]). 

In the current setting, we will simply set H.q = 7i. We presume that the equation 
<j2j) is given on the domain T>q®£ and that it satisfies conditions HHH] We also 
presume that T>o = span{?/v : £ £ where {ipe}e>o is an orthonormal basis of 

Ti. For simplicity, let us assume that N%j = S%f, this is not essential, see remark [7] 
below. 

Define PW to be the orthogonal projection onto span{-0o, . . . , ^fc}- We proceed 
to define the equations (P) by truncating the coefficients. Since the truncated 
operators will be bounded, we set V = H. and we letMf^if andifW in equation 
(HJ) be given by 

L[ k) = PW^PW, M[ k) = -Zp, = PWJTPW 

for all 1 < i < n. Note that condition Q] is thus satisfied, and as the coefficients are 
bounded condition [5] is as well. Conditions [T]-[3] are satisfied for the limit equation 
by assumption. 

Proposition 2 (Finite dimensional approximation). Under the above assumptions, 
the truncated equations (Qp converge to the exact equation (0): 

lim sup || - U?i(i\\ =0 for all ^ £H <g> T. 

fe^oo o<t<T 

Proof. By theorem[2l it suffices to show that for every a, j3 £ C™ and u £ span{^f }i>o, 
there exists a sequence {u^ k '} C H such that 
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We may simply take u^ k ' = P( k 'u. Since u is an element of the linear span of 
{ipe}e>o, there is a C G N such that = u and ^, a f)) u (k) = %{<x&) u for all 
k>C. □ 

Remark 7. The assumption that iVy = Sij is not necessary for the result to hold; 
however, care must be taken to approximate N t j in such a way that the Hudson- 
Parthasarathy relations of condition [2] remain satisfied (simple truncation as above 
typically does not achieve this purpose) . Similarly, it is not difficult to show that the 
result still holds if only span-j^ : £ G Z + } C T>q, provided that span{?/^ : I G Z + } 
is a core for «£?( a/3 ), a, f3 € C™. The chief technical advantage of this simple result 
compared to, e.g., an application of the results in |21j . is the very strong nature of 
the convergence. 

4. Proof of the main convergence theorem 

Before turning to the proof of the main convergence theorem, we provide a proof 
of the semigroup properties of T^ a ^ , lemma [T] 

Proof of Lemma\^ As Ut is a contraction cocycle, it follows easily that T^' is a 
strongly continuous contraction semigroup. As otl[o,t] 6 ® f° r every a G <C n and 
t > (recall that we assume that & contains all simple functions), we obtain from 
|(TJ for u, v G V 

(u <g e(al[ ,t\), U t v <g> e(/3/ [0 ,t])) = 

,a*Pt 



e 



u,v)+ I e- a ' Ps {u®e{aI [oA ),U s {^ (a0) v®e{l3I [oA )))ds 



where we have written for u E T>q 

(n n n 

ij—l 2—1 i— 1 

Therefore, we obtain for it, v G T>o using the chain rule 

(4) (u, T t {aP) v) = (tt, v) + / (u, Tj Q « ^ a ^v) ds, 



with %W>v defined as in equation ([3]). Using the fact that T>q is dense in 7io and 
that 

\{u,T^ ^ a ^v)\ < \\u\\ \\T^^f (a ^v\\ < \\u\\ \\^ a ^v\\, 
we find using dominated convergence that Q holds identically for all u E Ho, 
v G D , and i > 0. Moreover, as ||T t (Q/3) «5?( a/ %|| < ||_Sf ("^H is bounded, it follows 
that T^'^C^^v is Bochner integrable and thus evidently (see, e.g., [HI section 
1.2]) 



T [ t aP) v = v+ ^ a ^v ds Vt > 0, v G V . 



But then we obtain using the strong continuity of T t 



( a /3) 



It v — V 



t\o t 



which establishes the claim. The result for T t ' follows identically. □ 
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The remainder of this section is devoted to the proof of theorem [2] We first 
prove weak convergence of f7 t to Ut for every t < oo; this is the content of 
the following lemma. The weak convergence will subsequently be strengthened to 
strong convergence, and ultimately to strong convergence uniformly on compact 
time intervals. 

Lemma 8. If for every a, (3 € C™ and u <E T>o, there exist u^ k ' € Dom(Jzf ( fe ' Q ^)) so 
that 

then (ipi, Ut'ip2) fc ^°°> (ipi, Utip2) f or every ipi,tp2 G Ho ® J 7 an d t < oo. 

Proof. As Ut and are adapted, it suffices to restrict our attention to Ho <X> T t y 
Let W C Ho® T t \ be a total subset. Then for any ipi,ip2 € "Ho ® T t \ and e > 0, 
there exist ^4,^4 G spanW such that \\ipi — ip[\\ < e, \\1p2 — ^2 II < e > an d 

UmsupKVi,^-^)^)! 

< 2 ||Vi - HV2II +2 ||V>2 - ^11 Mil + limsup (C/ t (fc) - C/ t )^ 2 )l 

A;— »oo 

< 2 £ (||^|| + ll^ll) + 2e 2 +limsup|(^ 1 ,(C/ t (fc) - C/ t )^)l- 

— >oo 

It thus suffices to prove that the limit on the right vanishes for every 7/4, ip' 2 € spanW, 
i.e., it suffices to prove the result for ^1)^2 G W only. 

Now consider the total subset W = {w <£> e(f ) : u 6 2?o, / G ©'}> where & 
is the set of simple functions in L 2 ([0,i];C") (recall that by assumption & is 
admissible, so that U C X>o Let ipi £U, i 6 {1,2}. Then there exist Ui GT> , 

= t < ti < ■ ■ ■ < t m < t m+ i = t, and Oio, • • • , «im G C™ such that ipi = «i®e(/i) 
with fi{s) = cxij for s € [tj,ij+i[- It is not difficult to verify that, by virtue of the 
cocycle property, 

C/ t (fe V 2 ) - ||e(/i)|| ||e(/ 2 )|| (ui,r* aioa " o) T4^ 1 " a "J • ■•T^^\ 2 ), 
and similarly for (?/>i, Ut^-i)- In particular, the result follows as 

|(Vi,(^ (fc) - W 2 )l 

= l|e(A)|| ||e(/ 2 )|| |<ui, (T ( t ^ awa2a) ■ ■■T^ a2j) - T^ waw) ■ ■ ■ T ( t "\^ ] )u 2 ) \ 
< ||e(/i)|| ||e(/ 2 )|| |M| || (T^' awa2o) ■ ■■T t ^ t y a23) - T ^ loa2a) ■ ■■T^ 2i) )u 2 \\, 
which converges to zero as k — > 00 by the Trotter-Kato theorem. □ 

Proof of Theorem^ (a)^(b). Weak convergence of C/ t (fc) to U t was proved in the 

(k) 

preceding lemma. But note that as U{ are contractions and Ut is unitary, we can 
write 

WW-UMf = ((ui k) ~U t )^,(ui k) -U t m < 2Uf -2Re(([/f ) >, U&)). 

As (if), u\ Uftp) — > \\ip\\ 2 as fc — > 00 follows from the previous lemma, we immedi- 
ately obtain strong convergence. It thus remains to prove that strong convergence 
for every fixed time t can be strengthened to strong convergence uniformly on com- 
pact time intervals. To this end, we will appeal again to the Trotter-Kato theorem 
in a slightly different manner. 
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It is convenient to extend the Fock space to two-sided time, i.e., we will consider 
the ampliations of all our operators to the extended Fock space T — T s (L 2 (R; C n )) = 
■7-L<X> T , where T- = T is the negative time portion of the two-sided Fock space. We 
now define the two-sided shift § t : L 2 (R; C n ) -> L 2 (R; C") as 9 t f(s) = f(t + s), and 
by Qt '■ J~ T its second quantization. Note that Q t is a strongly continuous one- 
parameter unitary group, and that the cocycle property reads U t + S — U s Q* s U t Q s , 
etc., in terms of the two-sided shift. Now define on the two-sided Fock space the 
operators 

v t {k) = e t £/ t (fc) *, Vt = ® t u;. 

Then it is immediate from the cocycle property that and V t define strongly 
continuous contraction semigroups on Ji <£> T and 7io ® respectively, whose 
generators we will denote as JzfW and (see (TB] for a detailed study in the 
bounded case). Moreover, 

IIV-II 

as 0t is an isometry, where ip- ® V S .F_ ® W ® F = H. <g) .F. Hence, by the 
Trotter-Kato theorem, it suffices to show the following: for any ip G Dom(«5f), there 
exist V> (fc) 6 Dom(«5f( fe )) such that i/;( fe ) -> ^ and jSf( fc )^W -> jSfty as fe -> oo. Let 
us prove this assertion. Fix -0 G Dom(Jzf) and A > 0, and define 

poo 



ip w = / e" At F t w - JSfty) dt = R\ > (Xip - 
Jo 

where R^' is the resolvent of v} k ^ . We have already shown that || U^*ip — U^ip\\ 
as k — > oo for every fixed t > and ip G 7io ® -T 7 , so that evidently 

Therefore, we find by dominated convergence that 

poo 

^ (fc) J fe =! oo^ / e -*t Vt (\ip - &ip) dt = R x (\ip - Jfip) = ip, 



where we have used a standard result on resolvents section 2.1]. Similarly, we 
find that 

tpW G Dom(y( fc »), if( fe )^ (fc) = J£ {k) R {k) {Xip - JSfty) = JSfty + A(V> (fc) - V) 
for every fc G N by virtue of another standard result on resolvents [5J loc. cit.] . Thus 
evidently jSf WV (fc) -^V> 

as fc — > oo, and the proof is complete. □ 

Proof of Theorem^ (b)=>(a). Choose a,/3 G C n , v G 7Y , and let ip = v(g>e(pi[ Q . t ]). 
We obtain an estimate on || (T^ k ' a ^ — T^ a ^)v\\ through the following steps: 

||(7f ;Q/3) - T t {aP) )v\\ = sup |(u, (lf ;Q/3) - T ( t aP) )v)\ = 



||n||<l 

\(u®e(aI [oA ),(U^ -UtM\ \(iP\(ut k) -U t )iP}\ 

™£ e (M 2 +l/3| 2 )*/2 ",™ P ^ el/3! 2 */2 

IMI<i Wll<i 



-1/31^/2 || (f/ W _ W || < e -l/slVy||^|| a _2He«I7 t (fc ^,^)). 
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But by assumption ||(J7 t — P t *)?/>o|| —* as k — > oo for all tpo £ H.q J 7 , so that 
in particular (^),u[ k ^*Uti>) — > HV'II 2 as fe ^ oo. Wc thus obtain 

||(lf ;a/3) -T t ( ^>|| Va,/3 G C", v G Ho, t < oo. 

Now denote the resolvents R<£ ;Q/3) = (A - jSf(*;«0))-i. We fix A > 0, u G £>o, and 
define 

/* oo 

tt<*> = / e- At r t (fc;Q/3) (A-^ Q/3 )) U ^ = 4 fci Q «(A- JS?< a «)u. 
Jo 

Then it follows as in the proof of (a)^(b) that u< fc ) -► it and jSf(*>«0) u (*) jgf 
as fe — > oo (in particular, u( fe ) G Dom(Jz? w"/ 3 ))). This establishes the claim. 

It remains to show that if V is a core for all J£( k '< a P\ k £N, then we may choose 
{u^} C V. Indeed, let us fix a,/3 £ C™, u £ V , and construct the sequence 
as before. As V is a core for j£f we find that {(w, if (A;;Q/5) w) : t> € 2?} is dense 
in Graph(Jzf C '"' 3 )). Therefore, for every k, we can find a vector w^ k ' £ V such that 

|| w (fe) _ u (*)|| < \\£>(k;a0) w (k) _ ^(fe;a/3) n (fc)|| < fe -l 

Then w< fe ) -> u and £>(k-, a p) w (k) _> _^ (<*/3) U) and thc proo f is com plete. □ 

5. Proof of the singular perturbation theorem 

Before devoting ourselves to the proof of theorem [3l let us show that the coeffi- 
cients in assumption [3] satisfy the unitarity relations in condition [2] 

Proof of Lemma® We must verify the following relations: 

n n n n 

i=l j=l j=l j=l 

on T>q. The problem is that products such as Lih\ will contain terms such as 
GiPoG\, which can not be related directly to the unitarity conditions for K^ k \ 
etc. The additional conditions in assumption [3] are chosen precisely to alleviate this 
problem, as they allow us to cancel the Pq terms sandwiched inside the products. 
For example, as P {G t - AYF t )Pi = 0, 

n n 

Y^UL\ = - AY Fl )(G\ - F/ytAt)^. 

i=l i=l 

Note that, by applying condition Q] and assumption Q] to and K^ k \ 
k 2 {Y + yt) + k(A + At) + B + Pt = + = 

n n 

i=l i=l 
and as this must hold for all fc we obtain 

n n n 

Y+Y^ = ~Y F i F h A+A^ = ~Y( F i G i+ G i F i)> B+B^ = -Y G i G l 

i—1 i—1 i—1 
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Note, in particular, that as F- Pq = (assumption Efd)) 

n n 

(A + At)P = -J2 FiG\P , P (A + At) = - £ PoG,^. 



j=i 

Thus 



L,L\ = P \(A + At)Y"Ut + AY (A + At) — P — pt _ (y + yt)f t^t 



Pn 



Using that PoAPo = (assumption dfc)), it follows that this expression coincides 
with —K — K' as required. We now proceed to proving that the relation for Mj 
holds. Note that 

n n 

J2 N v L ] = E p oWu(S tj + FjYF^G] - f]y^A^)P , 

3=1 j,e=i 

where we have used assumption (J3J) as above. Using again assumption ^ and the 
above relations, we obtain through straightforward manipulation 

n 

Y,Vti + fIyFjKg] - F J tyt A t)p = (G t _ f }ya)p , 

3=1 

so that the claim follows directly. Next, we turn to the relations for iVy. We must 
establish 

n 

£ PoW ma (S aj + F^YFjWbj + FjY^F b )w} b P = 5 me 

j,a,b=l 

on l?0j where we have used assumption ([3]). Using assumption @ and the above 
relations as before, we obtain through straightforward manipulation 

n 

^(5 aj + FlYF,)^ + F^F b ) = <W 

3=1 

Hence the unitarity relation for Wij establishes the claim. It remains to show that 

n 

Poihi + F^F b )wj b W ja {S am + Ptyp m )p = S mt . 

j,a,b=l 

Using the unitarity relation for we find that we must show 

n 

Y Po{S a e + F}Y^F a )(S am + FlYF m )P = 5 mt . 

a=l 

But this was already established above, and we are done. □ 

We now turn to the proof of theorem [3J We will apply the general convergence 
result, theorem [2 which leaves the problem of finding the appropriate sequences of 
vectors u^ k \ To this end, we employ a simple but cunning trick due to Kurtz [19) . 

Proof of Theorem^ By theorem[2j it suffices to show that for every a, (3 € C n and 
u e Vq, there exists a sequence {u^} C T> such that 
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Let v € T>, and note that in the current setting ((3]) reads 
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i I fl|2 



£ a*W l3 P 3 + Gift - Y, a tWijGl +B- LL^L. 



:= (fc 2 y + fc^ Q/3 > +5^)?;. 
We will seek of the form i^ fc ) = u + k~ l ui + k~ 2 U2 with ui, «2 S I?. Note that 
g(k;ocf}) u {k) = k 2 Yu + k{Y Ul + A {aP) u) + Yu 2 + A {aP) Ul + B {aP) u + o(l). 



Evidently the proof is complete if we can choose U\, U2 in such a way that ^ a P) 

u + o(l) for all a, ft £ C n and u € T>q. The structural requirements (assump- 
tion [2]) are chosen precisely so that this is the case. First, note that Yu — by 
assumption Hfb), as u € T>q. This takes care of the quadratic term. Next, the linear 
term (oc k) must clearly vanish. Hence it must be the case that 



spa,n{Au,WijFjU, FiU : u S X>o, i,j = 1, . . . , n} C {Yv : v e T>}. 



But note that WijFju — for all u £ Vq by assumption HJd), while assumptions 
M.d,e) imply that F, = P x F t and AP Q = P^APq. Thus F t v = YYF t v and Aw = 
YYAv for all u € £>o by assumption^c), which establishes the claim. In particular, 
we will choose Ui = —YA^'u, which cancels the linear term. It remains to choose 
U2 € T> so that 



Note that even if we did not know in advance the form of K, Li, Mi, Nij, we must 
require that ££^ a ^u 6 7io for every a £ Do in order that the limit equation leaves 
H.0 invariant. Hence, if the limit of the slow dynamics exists, we must be able to 
choose U2 so that 



P\Yu 2 = -Pi{B^ - A^YA^y 
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But P\Y = YYY = YP\ = Y by assumptions [2^b,c), so that evidently we must 
have 

jSf(«0 w _ p Q (s(«/3) _ A {af ^YA {ap) )v = 

I n / n n \ 

p ° { E a u E ^(^^- + <%) + E ^-^^ Pi 

\£=1 £=1 ) 

71 

< {Wn{G\ - FjYA) - AYWijFj] 



+ > ((/, .IV- - F, YA) S, -tB ay a 

i=l 



n 



E F ^ 



y 



E F ^ 



E «i 



y 



E a * w v F i 



for all v € 2?0j ct, /? € C™. It remains to verify two things: that this expression does 
indeed define, for every a, /3 € C™, the generator of the semigroup and that 

y is invertible on the range of P X {B^ - A^YA^)P . But the latter follows 
immediately from assumption [2jc), so it remains to deal with the former. 

In order for the above expression to define the generator of , at the very 
least the last two terms must vanish — after all, they are quadratic in /3j and a*, 
respectively, which is inconsistent with ([3]). However, this is immediate from as- 
sumption HJd) , as it implies that FjPq = PoFi — 0. Simplifying further, we find 
that for all v S T>q 



J? iap) v = Po<| E ^WuiFjYFj + - £ ^W tJ (G] - F}YA) 



■J2(Gi - AYFApi + B- AY A - 



Thus evidently the coefficients of assumption [3] emerge naturally from our approach 
(by lemma Q} , and the proof is complete by lemma [2] and conditions HHSl □ 
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